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Initial Helioseismic Observations by Hinode/SOT 



O 
O 

(N 



6 



> 

O 

00 

o 
o 



X 



Takashi Sekii^, Alexander G. Kosovichev^, Junwei Zhao^, Saku Tsuneta\ Hiromoto Shibahashi^, Thomas E. 
Berger"^, Kiyoshi Ichimoto\ Yukio Katsukawa\ Bruce W. Lites**, Shin'ichi Nagata^, Toshifumi Shimizu^, 
Richard A. Shine'', Yoshinori Suematsu^, Theodore D. Tarbell'', and Alan M. Title'' 
^ National Astronomical Observatory of Japan, Mitaka, Tokyo 181-8588 
sekii@solar.mtk.nao.ac.jp 
^ W. W. Hansen Experimental Physics Laboratory, Stanford University, Stanford, CA 94305, USA 
^Department of Astronomy, School of Science, University of Tokyo, Bunkyo-ku, Tokyo 113-0033 
*High Altitude Observatory, National Center for Atmospheric Research, P.O. Box 3000, Boulder, CO 80307, USA 
^ Kwasan and Hida Observatories, Kyoto University, Kuravashira, Kamitakara-cho, Takayama, Gifu, 506-1314 
^Institute of Space and Astronautical Science, Japan Aerospace Exploration Agency, 3-1-1 Yoshinodai, Sagamihara, Kanagawa 

229-8510 

''Lockheed Martin Solar and Astrophysics Laboratory, B/252, 3251 Hanover St., Palo Alto, CA 94304, U.S.A. 

(Received 2007 July 27; accepted ) 
Abstract 

Results from initial helioseismic observations by Solar Optical Telescope onboard Hinode are reported. 
It has been demonstrated that intensity oscillation data from Broadband Filter Imager can be used for 
various helioseismic analyses. The k — uj power spectra, as well as corresponding time-distance cross- 
correlation function that promises high-resolution time-distance analysis below 6-Mm travelling distance, 
were obtained for G-band and Call-H data. Subsurface supergranular patterns have been observed from 
our first time-distance analysis. The results show that the solar oscillation spectrum is extended to much 
higher frequencies and wavenumbers, and the time-distance diagram is extended to much shorter travel 
distances and times than they were observed before, thus revealing great potential for high-resolution 
helioseismic observations from Hinode. 
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1. Introduction 

Local helioseismology has been a great success in reveal- 
ing subsurface structure of the sun, particularly in detect- 
ing three-dimensional flow and sound-speed anomaly in 
active regions (see, e.g., Zhao, Kosovichev & Duvall 2001, 
Zhao & Kosovichev 2003). 

Observations for local helioseismology require high res- 
olution and wide coverage both in space and in time. The 
50-cm aperture Solar Optical Telescope (SOT, Tsuneta 
et al. 2007; Suematsu et al. 2007) onboard Hinode satel- 
lite (Kosugi et al. 2007) provides 0.2-arcsec spatial res- 
olution though with a limited field of view (328 arcsec 
xl64arcsec for Narrowband Filter Imager (NFI) and 
218 arcsec x 109 arcsec for Broadband Filter Imager (BFI) , 
nominally). Except during the satellite eclipse periods, 
Hinode can observe the sun continuously with cadence up 
to 15 s or even less, with the actual limit placed by teleme- 
try bandwidth and capacity of the data recorder, which 
are shared with other instruments as well as other SOT 
programs running at the same time frame. 

The high resolution of SOT particularly offers a great 
potential for local hclioscismology. One of the main sci- 
entific targets of Hinode is to investigate the generation, 
evolution and eventual dissipation of active regions on the 
sun, in attempt to observationally constrain theories on 
dynamo processes which are thought to be responsible for 



the solar activity cycle. Local helioseismology can con- 
tribute greatly in this area, not only via active region to- 
mography but also in measuring flows at larger scales such 
as meridional flow, which plays an important role in the 
current models of solar activity cycle (Dikpati & Oilman 
2006). Please note that even with the limited field of 
view of SOT, and hence the limited depth coverage for to- 
mographic study (see below), one can still measure near- 
surface components of large-scale flows in a small area 
flrst but then extend the measurement to a broader area 
simply by repeating the observation at different locations. 

It was planned to use Dopplergrams produced by NFI 
for helioseismic observations on Hinode (Sekii et al. 2001; 
Sekii 2004). Not only oscillation signals in Dopplergrams 
are known to exhibit higher signal-to-noise ratio than os- 
cillation signal in intensity, but also NFI has a larger field 
of view than BFI. A larger field of view translates into 
a better depth coverage for tomography, because waves 
that penetrate deeper into the sun also travels further in 
horizontal directions. Thus NFI Doppler observation was 
the preferred mode of helioseismic observation by Hinode. 

However, long series of Dopplergrams are still only be- 
ing tested at the time of writing. Instead, we examined 
BFI instensity oscillation data from H line (Call) and G 
band (CH). 

The main aim of the present paper is to report that 
helioseismic observations can indeed be made with BFI. 
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As an example, a result from our first attempt at time- 
distance analysis of subsurface flow is presented. 

2. Data analysis 

A quiet region around the disk center was observed from 
around 1600UT, 1 Jan 2007, to around 0400UT, 2 Jan 
2007, for about 12 hours with a 12-minute interruption 
in the middle of the run. Hinode was tracking the re- 
gion during the observation (except during the interrup- 
tion mentioned in the above) and in addition Correlation 
Tracker was used to maximize framc-to-framc correlation 
in the central region of the field of view (Shimizu et al. 
2007). No further tracking was done on ground, unlike in 
the case of a sunspot study (Nagashima et al. 2007). 

Both G-band images and Call-H images were acquired 
at approximately 1-min cadence, Call-H images following 
G-band images by 20 s. Sample images are shown in figure 
1. The CCD pixels were 2x2 summed onboard to reduce 
the data amount, to double the effective pixel size to 0.106 
arcsec. We applied another 2x2 summing on ground (the 
overall result was a 4 x 4 summing) for the sake of faster 
analyses, though in this case we no longer fully exploit 
the 0.2-arcsec resolution of the telescope. The benefit of 
using the full 2 x 2-summed data, or observing with the 
full capacity of SOT (no summing), will be investigated 
in future work. 

The cadence was not exactly regular, due to the fact 
that there are two clocks involved in controlling the 
equipment. Not only they are on different cycles (0.5 s 
for Mission Data Processor and 1.6 s by Polarization 
Modulator Unit), the clocks are slowly drifting away from 
each other. Although programming tricks to correct this 
behavior are being investigated and tested, they have not 
been applied to the current datasets. 

A simple test was run to check the effect of the slightly 
irregular cadence, by putting through a sine wave on this 
irregular temporal grid and then computing power spec- 
trum, pretending that the cadence was exactly 1 minute. 
As we expected, the effect was insignificant and therefore 
we carried out all the following analyses assuming that 
the cadence was exactly 1 minute. 

The projection effect was also neglected. Because the 
BFI field of view is so narrow it should not give rise to 
any significant inaccuracy around the disk center. 

We used only the first uninterrupted 6-hr of the data 
for basic diagnostics products in the next section. For 
the time-distance analysis, the full 12-hr data were used. 
Before starting the analyses that follow, first we applied 
running difference of images: each image was subtracted 
from the following image, thus obviating the need for cali- 
brating each image for dark, flat fleld and bad pixels. One 
might worry that employing the running difference might 
increase the noise level, if the cadence is too short com- 
pared to the time scale over which the signal changes, in 
which case the differences between two successive frames 
arc dominated by noise. This is not the case here, when 
the 1-min cadence is not too high for the oscillation signal 
that is dominated by 5-min oscillations. 
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Fig. 1. Sample images of the observed region in G band (top) 
and Call H line (bottom). Of the full field of view which is 
218 arcsec X 109 arcsec, only the central part of the left half 
(approximately 54 arcsec X 27 arcsec) is shown. 



3. Power spectra and time-distance diagram 

In this section we present the most basic products 
of hclioscismic observations: the so-called k-uj diagram 
(or l-v diagram) and the time-distance diagram. They 
are then compared with results obtained from Michelson 
Doppler Imager (MDI) onboard SOlar and Heliospheric 
Observatory (Scherrer at al. 1995). The MDI data we 
used were acquired in the high-resolution mode. Although 
they too cover the central region of the disk, they are not 
coeval with the SOT data. 

Figure 2 shows the k-uj diagrams obtained from G-band 
and Call-H data; the series of the differential intensity- 
field data (intensitygram) were Fourier transformed and 
converted to power spectra, which were then azimuthally 
integrated in wavenumbcr space. Although no spherical- 
harmonic decomposition was involved, degree I was as- 
signed to each component by multiplying its wavenumbcr 
by the solar radius. The power spectra were corrected for 
the factor, where cu is angular frequency, for we have 
taken running difference before Fourier transform. 

In spite of only 6 hours of observation, the f-mode ridge 
is seen up to I 4000. We also see the high-frequency in- 
terference peaks (Jefferies et al. 1988, Duvall et al. 1991, 
Kumar et al. 1990) up to the Nyquist frequency. They 
seem disturbed close to the Nyquist frequency because 
the ridges from both the super- and sub-Nyquist ranges 
are running into the other regime. Compared to the k — uj 
diagram obtained from a 512-min MDI high-resolution ve- 
locity dataset, also shown in figure 2, the better perfor- 
mance of SOT at high- wavenumbcr, where high-resolution 
matters most, is noticeable in how the f-mode ridge ex- 
tends to higher wavenumbers and, in the case of Call H 
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Fig. 2. Power spectra (k—uj diagrams) from the 6-hr G-band 
data (top) and Call-H data (middle), compared with power 
spectrum from a 512-min MDI high-resolution Doppler data 
(bottom). 



data, the pi ridge as well. However, as was expected, the 
SOT k — Lu diagrams show higher background-noise level 
compared to the Doppler measurement by MDI. 

It is apparent that in G-band data the convectivc 'noise' 
is stronger and for this reason, in the following analyses, 
we use only Call H data. More detailed comparison be- 
tween the two, in particular phase relations at various 
frequencies, would be a subject of future work. 

We then applied Fourier-Bcssel transform (see, e.g., 
Sekii & Shibahashi 2003) to the Call-H k-ui power spec- 
trum to obtain time-distance cross-correlation function 
(figure 3). The power spectrum uncorrected for the 
was used; this may be viewed as an application of a broad 
high-pass filter to the power spectrum, the intention being 
filtering out the convective noise. 

The structure down to less than 0.1 heliographical de- 
gree is seen with a high signal-to-noise ratio, which is very 
encouraging. Once again, we see the higher performance 
of SOT, compared to MDI, also shown in figure 3, at small 
spatial scales. The signal-to-noisc ratio, however, is higher 
for MDI at large distances and thus the both instruments 
complement each other. 

4. Subsurface flow 

We have analyzed the full 12-hr long scries by time- 
distance technique for local helioseismology. The gap 
between the two 6-hr stretches was filled with zeroes. 
Then for measuring travel-time difference between out- 
ward and inward wave components, each pixel was cross- 
correlated with average signals over annuli that cover, in 
radii, 4.00-6.62, 6.47-9.70, 9.24-12.47, 12.01-15.86, 15.09- 
18.94, 18.48-22.33, and 21.56-25.41 (all in Mm). One such 
example as well as an MDI counterpart (in high-resolution 
mode) is shown in figure 4. The annuli were further di- 
vided into 4 sectors representing north, south, east and 
west, and opposing pairs of sectors were similarly cross- 
correlated to yield northward-southward and westward- 
eastward travel-time differences. 

All the travel-time differences were then inverted by us- 
ing ray-approximation kernel for acoustic wave propagat- 
ing through the solar interior. Figure 4 shows outward- 
inward travel-time difference map, as a proxy divergence 
of the flow field. Inversions for horizontal flow in the 
depth ranges of 0-1 Mm, 1-2 Mm, 2-3 Mm and 3-4 Mm 
are shown in figure 5. Inversions were done for the en- 
tire area in figure 4, but to avoid overcrowding the figure 
with arrows, in figure 5 only the central quarter is shown. 
Horizontal flow patterns are largely consistent with the 
outward-inward travel-time difference, as was seen in re- 
alistic numerical simulation by Zhao et al. (2007). This 
indicates that the travel-time difference is a good proxy 
of flow divergence, though to acquire additional details we 
do need inversion procedure. We see that supergranula- 
tion patterns are coherent vertically within a spatial scale 
of a few Mm. 
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Fig. 3. Time-distance diagrams from 6-hr Ca II-H data (top) 
and from 512-min MDI high-resolution Doppler data (bot- 
tom). The datasets are not co-eval. 
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Fig. 4. An outward-inward travel time difference maps from 
512-min MDI high-resolution Doppler data (top) and from 
Call-H data (bottom), indicating flow divergence due to su- 
pergranulation. The inner radius of the annulus used was 
13.86 Mm and the outer radius 15.86 Mm. Dark patches are 
where inward travel time is longer than outward travel time 
i.e. where the flow is generally diverging, whereas brighter 
patches are where the flow is generally converging. Since the 
two datasets are not coeval, individual features are not to be 
compared directly. 

5. Discussion 

As we have demonstrated, SOT observation of intensity 
oscillations can be used for high-resolution hclioscismic 
diagnosis of the sun, in spite of its noise level higher than 
what can be expected from Doppler measurement. This is 
due to the fact that solar noise is higher in intensity than 
in Doppler velocity. 

From figure 2 it is apparent that SOT offers good op- 
portunity to study high-wavenumber waves, particularly f 
modes. Ishikawa et al. (2007) has found ubiquitous small- 
scale horizontal magnetic field from Hinode/SOT observa- 
tion of the photosphere. In addition to turbulence, pres- 
ence of such magnetic fluxes would affect propagation of 
f-mode waves and this would be an important future work. 

In the time-distance cross-correlation (figure 3), the 
structure below 0.5 degree, which is about the distance 
of 6 Mm, has been observed much clearer than before. 
Sekii (2004) used G-band data obtained from 50-cm La 
Palma SVST to demonstrate somewhat similar enhance- 
ment in spatial resolution but the noise level was much 
higher. It was speculated that the noise came mainly from 
atmospheric seeing and that the spaceborne SOT would 
perform better. Our current results show that SOT indeed 
does. Figure 4 demonstrates that study of fine structures 
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in supcrgranulation patterns benefits from tlie liigh reso- 
lution of SOT; the bright and dark patclies have a typi- 
cal spatial scale of 20 Mm but a transition from bright to 
dark takes place over the scale of less than 10 Mm, which 
MDI have difficulty in resolving. With SOT results we see 
smoother transitions. 

To exploit this small-scale information would require 
more tuning and testing of the inversion, particularly 
the setup of annuli. The trade-off between resolution 
and error need to be examined carefully. Once done, 
we will certainly achieve higher horizontal resolution in 
time-distance inversions. A similar improvement in per- 
formance is expected in resolving, vertically, layers im- 
mediately below the surface, since SOT can use higher- 
wavenumber and hence shallower-penetrating waves. This 
applies also to sound-speed inversion, which has not been 
attempted yet. 

Time-distance analysis in active regions, be it for flow 
or for thermal (and even magnetic) structure, would need 
further care particularly if we wish to go a step beyond 
phenomenological inversions; wave propagation in magne- 
tized atmosphere has to be somehow taken into account. 
For the initial analysis of sunspot oscillation data, see 
Nagashima et al. (2007). 

Hinode is a Japanese mission developed and launched 
by ISAS/JAXA, collaborating with NAOJ as a domestic 
partner, NASA and STFC (UK) as international partners. 
Scientific operation of the Hinode mission is conducted 
by the Hinode science team organized at ISAS/JAXA. 
This team mainly consists of scientists from institutes in 
the partner countries. Support for the post-launch op- 
eration is provided by JAXA and NAOJ (Japan), STFC 
(U.K.), NASA, ESA, and NSC (Norway). This work was 
carried out at the NAOJ Hinode Science Center, which 
is supported by the Grant-in-Aid for Creative Scientific 
Research "The Basic Study of Space Weather Prediction" 
from MEXT, Japan (Head Investigator: K. Shibata), gen- 
erous donations from Sun Microsystems, and NAOJ inter- 
nal funding. 
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Fig. 5. Subsurface flow maps (indicated by arrows) obtained 
by time-distance analysis of 12-hr Call-H data. Depth ranges 
are - 1 Mm, 1-2 Mm, 2-3 Mm and 3-4 Mm (top to bot- 
tom). The longest arrow in the panels correspond to 0.51km 
s~^, 0.52km s~^, 0.30km s~^ and The field of view corre- 
sponds to the central quarter of the field of view in figure 
4. 
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